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RESEARCH  ACCOMPLISHMENT  IN  THE  FIRST  YEAR 


The  research  program  was  planned  in  four  stages  over  two  calendar  years 
starting  September  28,  1979.  The  first  and  second  stages  would  encompass  both 
the  investigation  of  free  surface  properties  of  the  semiconductors  in  an  altra- 
high  vacuum  (lO"^”  torr)  system  (UHV)  equipped  with  LEED,  Auger  electron  spectro¬ 
scopy  (AES)  and  ion-sputtering  capabilties,  and  the  design  and  assembly  of  the 
satellite  chamber  for  device  fabrication.  These  two  stages  would  be  performed 
at  the  same  time  without  interference  with  each  other.  The  third  stage  of  the 
work  would  encompass  the  study  of  the  waveguide  surface  structure,  formulation  of 
surface-scattering  theory  and  characterization  of  scattering  loss  parameters  in 
terras  of  the  surface  geometry  through  LEED  and  SEM  analysis  in  the  same  ultrahigh 
vacuum  system.  The  last  stage  of  the  work  would  encompass  the  fabrication  of  the 
oxide-free  electro-optical  devices,  the  measuring  of  the  optical  propagation 
losses  and  the  study  of  the  optical  characteristics  of  such  waveguides. 

The  first  and  the  second  stages  of  the  planned  research  were  completed. 
Currently  we  are  conducting  the  third  stage  of  the  research  plan.  The  satellite 
chamber  for  device  fabrication,  as  shown  in  Fig.  1,  was  designed  and  assembled 
to  the  main  system.  A  schematic  diagram  of  the  main  system  is  shown  in  Fig.  2. 

The  satellite  chamber  is  indicated  as  subsystem  2  in  Fig.  2.  Subsystem  1  in  Fig.  2 
is  a  plasma  reactor  designed  for  oxidation  research  of  GaAs  and  is  in  no  inter¬ 
ference  with  this  research  program.  A  photograph  of  the  completed  system  is  shown 
in  Fig.  3. 

The  free  surface  properties  of  GaAs  (110),  GaP  (110)  have  been  studied  by  LEED 
and  AES.  Results  of  the  studies  have  been  presented  in  two  different  scientific 
conferences.  Abstracts,  which  summarized  the  results  of  the  research  works,  are 
given  in  Appendices  A  and  B.  Abstract  of  the  first  paper  entitled  "Chain  Method 
of  LEED/MEED  Intensity  Calculation  for  Diatomic  Surface"  is  given  in  Appendix  A 
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SPECIMEN  PLUG 


Fig.  1 .  Top-view  schematic  of  the  satellite  system  attached  to  the 
main  chamber  (left).  A:  sorption  forepump;  B:  Vacion  UHV  pump; 

C;  bellows-coupled  flange;  D:  Isolation  valve;  E:  six-armed  cross 
F:  ionization  gauge;  G:  transfer  probe  magnet  actuator.  In  the 
schematic  the  transfer  probe  is  shown  completely  retracted,  with  the 
specimen  plug  located  at  the  center  of  the  six-armed  cross. 
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1^1  g.  2.  Schematic  diagram  of  the  UHV  system  equipped  with  LEED,  AES,  mass 

spectrometer,  ion  sputtering  gun  and  ellipsometer  for  surface  analysis 
Subsystem  1  is  the  plasma  reactor  which  is  currently  under  constructio 
Subsystem  2  is  the  ‘CS  device  fabrication  chamber.  constructio 
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LEED/AES  system  for  Investigation  of  surface  ami  interface  properties  of  semiconductors 
and  ceramic  materials.  Tl>e  capital  equipment  is  acquired  througi»  Uie  USE  equipment 


and  the  second  paper  entitled  "Atomic  Structure  of  GaAs  (110)  Face"  is  given  in 
Appendix  B.  Both  papers  are  being  prepared  for  publication. 


The  third  stage  of  the  research  work  is  currently  ongoing.  The  theoretical 
formulation  of  surface  scattering  theory  and  characterization  of  the  scattering 
loss  parameters  are  currently  being  developed.  The  problem  of  a  wave  propagating 
along  z  direction  in  a  rectangular  dielectric  waveguide  with  dimensions  a  and  b  as 
shown  in  Fig.  4  was  formulated  and  solved.  The  wave  functions  in  the  waveguide 
were  solved  using  predominantly  polarized  approximation.*  The  propagation 
attenuation  constants,  e.g.  the  loss  parameters,  were  calculated.  Details  of  the 
theoretical  formulation  and  calculation  is  given  in  Appendix  C. 

Citations  of  research  results,  which  are  presented  in  conferences  or  sub¬ 
mitted  for  publication,  are  listed  in  Appendix  D. 


i 


*D.  Marcuse,  Bell  Syst.  Tech.  J.,  3187  (1969). 


Fig.  Cross  section  of  a  rectangular  dielectric  waveguide 
with  dimensions  a  and  b.  n] ,  ri2>  n3.  n4  and  nc  are 
the  indixes  of  refraction  of  regions  1  to  5,  respectively. 
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APPENDIX  A 


The  following  paper  was  presented  in  the  Conference  on  Determination  of 
Surface  Structure  by  LEED,  ISM  Thomas  J.  Watson  Research  Center,  Yorktown  Heights, 
N.Y.,  June  19-20,  1980 

CHAIN  METHOD  OF  LEED/MEED  INTENSITY  CALCULATION  FOR  DIATOMIC  SURFACES 

In  low  energy  electron  diffraction  (LEED),  crystal  surface  is  probed  by  bom¬ 
barding  electrons  of  energy  E  £  200  eV  at  normal  incidence  to  the  surface  and  analy¬ 
zing  the  intensity  of  those  elastically  scattered  electrons  which  are  back  reflected 
from  the  surface,  whereas  in  medium  energy  electron  diffraction  (MEED)  due  to  an  in¬ 
crease  in  the  energy  of  incident  electrons  (200  eV  <  E  £  5  keV)  the  surface  sensiti¬ 
vity  of  the  technique  is  maintained  by  taking  an  oblique  angle  of  incidence.  Inter¬ 
pretation  of  LEED/MEED  intensity  spectra  requires  detailed  calculations  of  the 
diffracted  intensities  for  a  series  of  trial  models.  For  sucn  a  technique  to  work, 
a  rapia  and  accurate  method  of  calculation  is  needed. 

Layer-KKR  method  is  conventionally  employed  in  the  calculation  of  LEED  intensity 
spectra.  In  this  formalism,  the  crystal  is  divided  into  a  number  of  layers  parallel 
to  the  crystal  surface  and  scattering  calculations  are  split  into  two  parts,  intra¬ 
layer  and  interlayer.  In  the  intralayer  scattering  calculation,  due  to  the  assumed 
spherical  symmetry  of  the  atomic  potential,  an  angular  momentum  representation  is 
used.  This  involves  large  matrices  and  their  inverses  and  becomes  cumbersome  with 
increasing  energy  and  more  complex  surfaces.  In  such  a  situation  we  propose  to  use 
the  chain  method  of  intralayer  multiple  scattering  calculation. 

In  the  chain  method  the  two  dimensional  intralayer  multiple  scattering  calcu¬ 
lations  are  further  divided  into  two  one  dimensional  steps,  scattering  within  a 
chain  of  atoms  and  between  the  chains.  At  each  stage  of  calculation  the  electron 
wave  function  is  represented  in  terms  of  an  appropriate  set  of  basis  functions.  The 
scattering  by  an  atom,  a  chain  of  atoms  and  a  layer  of  chains  is  expressed  in  terms 
of  spherical,  cylindrical  and  plane  wave  representation  respectively.  As  we  move  from 
one  stage  to  another,  transformation  from  one  basis  set  to  another  is  carried  out. 

Due  to  its  one  dimensional  lattice  summations  the  chain  scattering  formalism  has 
many  computational  advantages  over  the  layer-KKR  method  both  for  normial  incidence 
LEED  and  off  normal  incidence  MEED. 


APPENDIX  B 


The  following  paper  is  submitted  for  presentation  in  the  Eighth  Annual  Con¬ 
ference  on  the  Physics  of  Compound  Semiconductor  Interfaces,  Williamsburg,  Virginia, 
January  27-29 ,  1981 . 


ATOMIC  STRUCTURE  OF  GaP  (110)  FACE 
ABSTRACT 

Low  energy  electron  diffraction  (LEED)  intensities  have  been  measured  for  the 
(110)  face  of  GaP  and  analyzed  using  a  dynamical  multiple-scattering  model  of  the 
diffraction  process.  The  intralayer  multiple  scattering  is  treated  exactly,  while 
for  the  interlayer  multiple  scattering,  the  renormal ized- forward-scattering  method 
is  used.  Comparison  of  the  calculated  and  observed  LEED  intensities  suggests  that 
both  the  Ga  and  the  P  atoms  on  the  (110)  face  may  exhibit  a  contracted  outermost 
layer  spacing.  The  surface  l^y^C  ''s  compressed  by  about  5%  such  that  the  top  layer 
spacing  is  reduced  by  0.1  ±  0.02A.  The  rippled  geometry  of  surface  reconstruction 
is  not  clearly  observed.  This  indicates  that  the  GaP  (110)  surface  atomic  struc¬ 
ture  is  different  to  that  of  GaAs  (110).  By  comparing  the  structures  and  properties 
of  the  GaP  (110)  face  and  the  GaAs  (110)  face,  it  is  concluded  that  the  GaP  (110) 
face  is  relatively  unstable  and  reactive. 


APPENDIX  C 


The  problem  of  a  wave  propagating  along  z-direction  is  a  rectangular  dielectric 
waveguide  with  dimension  a  and  b  as  shown  in  Figure  4  is  formulated  and  solved.  An 
exact  analytical  treatment  of  rectangular  waveguide  is  practically  impossible. 
Therefore  the  approximated  analytical  approach  developed  by  Marcatili  (a)  is  followed 
in  solving  the  problem. 

Assuming  a  traveling  wave  propagating  in  z-direction  as  shown  in  Fig.  4, 

Maxwell  equations 


7XH  =  II 


can  be  solved  to  give 

3H 


0  3y 


H  = 


E,  = 


H  * 


E.  * 


3E  3H 

[-i/(n2k^'-3^)]C8  ^ 


[-i/(n^k^.3^)]C3;i-a.V^^] 


3H  3H 


(-JL _ X) 


z  E^n^ioj  ^9x  3y 


3E  3E 

0 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


where  n  is  the  index  of  refraction  of  the  dielectric;  k  is  the  wavevector  in 
free  space;  oi  is  the  angular  frequency  of  the  wave  and  3  is  the  propagation 
constant  in  z-direction. 

Substituting  (2)  and  (4)  into  (5),  and  (1)  and  (3)  into  (6),  we  have 


3^E  3^E 

=0  . (7) 
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52l|  :i2u 


Two  different  types  of  modes  can  be  supported  in  the  waveguide^  ,  e.g., 


Epq  :  polarized  predominately  in  x-direction  and 


:  polarized  predominately  in  y-di recti on, 

HH 

where  p  and  q  are  positive  integers  indicating  the  modes  of  the  wave  propagating 
in  the  waveguide. 

eL  mode 
pq 

In  rogion  1  with  refraction  index  n^  as  shown  in  Fig.  4,  the  waveequation  (7) 
can  be  solved  to  give 


E2(x,y)  =  A  Cos  {x+5)  Cos  (y+n) 


with  =  0,  where  C  and  n  are  the  phase  factors  of  E^. 

Other  components  of  the  wavefunction  can  be  obtained  by  substituting  {9a) 
into  eqs.  (1 ) ,  (3)  and  (4) : 


H  =  -  Sin  k„  (x  +  E)  Sin  k^,  (c  +  n) 


*^o  X 


(n?k^-k^  ) 

Ex  =  ■'A  - —  Sin  k^  (x  +  E)  Cos  k^  (y  +  n) 


Ey  =  — Cos  k^  (x  +  5)  Sin  k^  (y  +  n) 


Substituting  (9a)  into  (7)  we  have 
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For  small  incident  angle,  S»(k^  +  kM,  eq.  (10)  can  be  approximated  to  be 

X  y 

3  =  n^. k  Cos0  =  n^k  for  i<S°. 

The  wavefunctions  in  regions  2,  3,  4  and  5  with  indices  of  refraction  n^, 
n^,  n^  and  n^,  respectively,  as  shown  in  Fig.  4,  are  solved.  When  boundary  con¬ 
ditions  are  matched  between  region  1  and  the  neighboring  regions,  i.e.,  regions  2, 
3,  4  and  5,  we  obtain  the  four  transdental  equations; 


tan  k^a  =  n^k^  (n^yg  +  n^Y3)/(n|n|  kj  -  n:[Y3Y5) 

tan  k|  =  (k^/Yj) . 

tan  kyb  =  ky  (y2  +  Y4)/(k‘  -  Y2Y4)  . 

tan  k^n  =  -Y^/k^  . 

where  Y2  =  C(n|  -  n^)k^  - 

V3  =  [(nf  -  n|)k^  - 


(lla) 

(lib) 

(12a) 

(12b) 


Y5  =  [(n?  -  n5^)k^  - 

A  computer  program  is  currently  being  developed  using  Newton's  method  to  solve 

equations  (11)  and  (12)  to  obtain  k  ,  k  ,  6,  and  the  power  attenuation  constant  a. 

X  y 

REFERENCES: 

(a)  Marcatili,  E.  A.  J.,  Bell  Syst.,  Tech.  J.  48,  2071-2102  (1969). 

(b)  0.  Marcuse,  "Theory  of  Dielectric  Optical  Waveguide",  Academic  Press,  New 
York,  1974. 

(c)  I.  P.  Kaminow,  W.  L.  Mammel ,  and  H.  P.  Weber,  Applied  Optics,  Vol  13,  No.  2/ 
Feb.  1974. 
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Publication  citations  and  presentations  supported  or  partially  supported  by 
this  research  grant: 

1.  "Decomposition  of  Aluminum  oxide  by  Electron  Bombardment."  B.  W.  Lee  and 
J.  M.  Kuo,  B.  Am.  Phys.  S.,  25  (3),  238,  1980. 

2.  "Study  of  MIS  Polycrystalline  Silicon  Solar  Cell  Using  Auger  Electron 
Spectroscopy,"  J.  M.  Kuo,  B.  W.  Lee,  B.  Lalevic  and  W.  A.  Anderson,  B.  Am.  Phys.  S., 
^  (3),  409,  1980. 

3.  "Chain  Method  of  LEED/MEED  Intensity  Calculation  for  Diatomic  Surfaces," 

N.  Masud,  C.  G.  Kinni burgh,  0.  J.  Titterington.  Presented  in  Conference  on  Deter¬ 
mination  of  Surface  Structure  by  LEEO,  IBM,  T.  J.  Watson  Research  Center, 

Yorktown  Mights,  NY,  June,  1980. 

4.  "Study  of  MIS  Silicon  Cell  by  ESCA  and  AES,"  Y.  S.  Wang,  H.  J.  Yu,  C.  C.  Hsu 
B.  W.  Lee  and  W.  A.  Anderson,  presented  in  the  27th  National  Symposium  of  American 
Vacuum  Society,  Detroit,  October,  1980.  Submitted  to  J.  Vac.  Sci.  Technol  for 
publication. 

5.  "Stability  of  MIS  Silicon  Solar  Cell,"  B.  W.  Lee,  J.  M.  Kuo,  B.  Lalevic  and 
W.  A.  Anderson,  submitted  to  J.  Vac.  Sci.  Technol.  for  publication. 

6.  "Atomic  Structure  of  GaP(llO)  Face,"  B.  W.  Lee,  R.  K.  Ni  and  N.  Masud, 
submitted  for  presentation  in  the  Eighth  Annual  Conference  on  the  Physics  of 
Compound  Semiconductor  Interfaces,  Williamsburg,  Virginia,  Jan.  1981. 
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A  NONLINEAR  MAXIMUM  ENTROPY  METHOD  FOR  SPECTRAL  ESTirvATION 

Summary 

An  intensive  research  study  is  made  of  the  nonlinear  maximum 
entropy  spectral  analysis  method  proposed  by  P.  F.  Fougere  of 
AFGL.  The  research  not  only  provides  a  much  better  understanaing 
of  the  properties  of  the  method  including  spectral  resolution, 
convergence,  etc.  but  also  shows  that  the  method  can  be  perform'.'! 
effectively  with  the  PDP  11/45  minicomputer.  The  successful 
implementation  of  tne  method  is  described  in  the  program  listings 
in  Appendix  I.  Extensive  computer  results  are  presented  on 
various  data.  These  results  clearly  confirm  that  the  nonlinear 
method  is  superior  to  the  Burg's  maximum  entropy  spectral  analysis. 

Preliminary  results  on  the  multichannel  (multivariate)  maximum, 
entropy  spectral  analysis,  the  two-dimensional  maximum  entropy 
spectral  analysis,  and  computer  graphics  for  the  spectral  display 
are  also  presented.  An  extensive  bibliography  of  the  miaximum, 
entropy  spectral  analysis  is  given  in  Appendix  III. 
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I.  Review  of  Research  Progress 

This  research  is  concerned  with  the  nonlinear  maximu;n  entropy- 
spectral  analysis  (MESA)  method  proposed  by  Dr.  Paul  Fougere  of 
AFGL.  As  verified  in  this  research  work  performed  at  the 
PDF  11/45  minicomputer,  the  method  not  only  provides  a  much  better 
spectral  resolution  than  the  Burg’s  method  but  also  removes  the 
line-splitting  and  frequency  shifting  phenomena  for  sinusoidal 
signals,  as  experienced  in  the  Burg’s  method.  By  using  tr.e  double 
precision,  the  minicomputer  results  are  reasonably  close  to  those 
obtained  at  the  GDC  6600  computer  at  AFGL. 

The  final  computer  program  developed  for  the  nonlinear 
complex  sigrial  maximum  entropy  spectral  analysis  is  chowr.  in 
Appendix  lA .  The  program  follows  the  mathematical  development  of 
Fougere  [1]  but  is  quite  different  from  the  original  computer  pro¬ 
gram  developed  by  fr.  Fougere.  Appendix  IB  is  the  computer 
program  for  Burg's  complex  signal  maximum  entropy  spectral  analysis. 
By  using  the  computer  programs  and  the  two-channel  radar  aata  as 
shown  in  Fig.  1  (see  also [”2]  ),  the  spectrum  of  the  nonlinear 
method  is  shown  in  Fig.  2a  (linear  plot)  and  Fig.  2b  (logarithmic 
plot)  for  10  filter  weights.  The  Burg's  result  is  shown  in  Fig.  3a- 
(linear  plot)  and  Fig.  Jb  (logarithmic  plot),  also  for  10  filter 
weights.  More  detailed  tabulation  of  the  major  frequency  compo¬ 
nents  outside  the  clutter  bandwidth  (-0.l67fs,  +0.l67fs)  is  given 
in  Fig.  4.  The  nonlinear  method  which  is  very  close  to  the  true 
answer  clearly  is  much  better  than  the  Burg's  method.  30  iterations 
are  used  in  the  nonlinear  method  which  appears  to  be  an  optimum 
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number.  The  optimum  filter  weight  is  around  10  or  11  as 
determined  by  Fig.  5  which  shows  the  linear  prediction  error  (FiZ) 
according  to  the  Akaike  criterion.  The  solid  curve  is  the  lower 
bound  and  the  dashed  curve  is  the  upper  bound  as  tabulated  in 
Fig.  6.  The  upper  bound  was  originally  proposed  in  the  research 
proposal  [2]  and  it  seems  to  be  better  than  the  lower  hound  due  to 
Akaike . 

Major  documentations  already  made  which  describe  the  research 
progress  are  as  follows! 

1.  C.H.  Chen,  J.  Chen,  and  C.  Yen,  "A  minicomputer  implementation 
of  Fougere's  maximum  entropy  spectral  analysis  method,"  Techni¬ 
cal  report  prepared  for  the  Mini-grant,  August  20  1980.  This 
report  has  detailed  results  of  comparison  between  Burg's  and 
the  nonlinear  methods  for  sinewave  and  sunspot  data. 

2.  C.  H.  Chen,  "Spectral  resolution  of  Fougere's  maximum  entropy 
spectral  analysis,"  to  be  published  in  the  Proceedings  of  IEEE, 
June  1981.  This  journal  paper  based  on  the  work  performed 
under  the  Mini-grant  provides  a  good  comparison  between  the 
nonlinear  and  Burg's  methods,  and  the  other  method  for  complex 
sinusoids.  The  Cramer-Rao  bound  is  used  as  a  reference. 

3.  C.  H.  Chen  and  C.  Yen,  "Note  on  computer  graphics  for  maximum 
entropy  spectral  analysis,"  Technical  report  prepared  for  the 
Mini-grant,  March  23»  1981.  This  report  provides  a  three- 
dimensional  spectral  display  of  sinewaves  for  both  nonlinear 
and  Burg's  methods. 
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A  number  of  important  results  are  included  in  the  new  research 
proposal  submitted  to  AFOSR  in  December  1980.  Appendix  III 
provides  an  extensive  list  of  references  on  the  maximum  entropy 
spectral  analysis.  The  following  sections  describe  some  new 
research  areas  with  preliminary  results. 

II.  Multichannel  (multivariate)  Maximum  Entropy  Spectral  Analysis 

A  mathematical  presentation  of  this  topic  is  given  in  Appendix 

II.  Several  computer  programs  for  multichannel  maximum  entropy 
spectral  analysis  were  provided  by  Dr.  Fougere.  The  following 
results  are  based  on  the  time  series  of  sunspot  numbers,  northern 
light  activity,  and  earthquake  activity  by  using  the  third  multi¬ 
channel  program.  The  data  are  tabulated  in  [3]- 

Fig.  7a  is  the  first  channel  (sunspot  number)  auto-spectrum 
with  linear  (left)  and  logarithmic  (right)  scales,  and  l6  lags. 

Fig.  7b  is  the  second  channel  (northern  light  activity)  aumo- 
spectrum  with  linear(left)  and  logarithmic  (right)  scales,  and 
16  lags. 

Fig.  7c  is  the  third  channel  (eartnquake  activity)  auto¬ 
spectrum  v/ith  linear  (left)  and  logarithmic  (right)  scales  and 
l6  lags. 

Fig.  7d  is  the  cross-spectrum  between  channels  1  and  1  with 
real  part  (left)  and  imaginary  part  (right),  and  l6  lags. 

Fig.  7e  is  the  cross-spectrum  between  channels  2  arid  3  v/ith 
real  part  (left)  and  imaginary  part  (righ  ),  and  16  lags. 

Fig.  71  is  the  cross-spectrum  between  channels  1  and  3  v/ith 
real  part  (left)  and  imaginary  part  (right),  and  l6  lags. 

Dy  way  of  verification,  it  is  interesting  to  note  that  the 
spectral  peaK  for  the  sunspot  numbers  is  determined  accurately. 

III.  Two-Dimensional  Maximum  Entropy  Spectral  Analysis 

he  consider  a  very  simple  separable  case  in  the  tv/s-dimensional 
spectial  analysis.  The  signal  considered  is  sin(2Fx)  sin(2?ry). 

In  tnis  case  the  power  spectrum  is  the  product  of  the  power  spectra 
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of  sin(2T[x)  and  sin(2Try).  In  each  spatial  dimension,  the  power 
spectrum  can  be  determined  by  using  the  Fourier,  nonlinear  and 
Burg's  methods.  The  two-dimensional  spectra  are  shown  in 
Figs.  8,  9  and  10  respectively,  based  on  the  Fourier,  nonlinear, 
and  Burg's  methods.  The  nonlinear  method  clearly  is  much  better. 
Extension  of  the  above  procedure  to  a  more  general  two-dimensional 
spectral  analysis  is  not  possible.  Although  some  two-dimensional 
maximum  entropy  spectral  analysis  work  has  been  reported  (see 
Appendix  III),  the  success  is  very  limited.  Further  research  is 
much  needed. 

IV.  Conclusions  and  Recommendations 

The  nonlinear  maximum  entropy  spectral  estimation  method 
proposed  by  P.  F.  Fougere  has  provided  superior  spectral  estima¬ 
tion  over  the  Burg's  method  in  a  number  of  data  considered.  The 
computational  requirement  of  the  nonlinear  method  is,  however, 
significantly  higher.  Typical  number  of  iterations  needed  is  20 
to  30.  The  use  of  minicomputer  has  not  created  significant 
computational  problem  as  predicted.  It  is  our  strong  belief  that 
the  nonlinear  method  will  become  very  popular  in  high  resolution 
spectral  analysis  for  a  wide  range  of  applications  in  geophysics, 
sonar,  radar  areas,  etc. 

Recommended  further  studies  include  the  multichannel 
(multivariate)  maximum  entropy  spectral  analysis,  the  two- 
dimensional  nonlinear  maximum  entropy  spectral  analysis,  computer 
graphics  for  the  spectral  display,  signal  decomposition,  and 


signal  prediction  and  extrapolation. 
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No.  of  weight 

Log.  of 

n+Tm+i)  _ 

Log.  of 

N+2(M+1 ) 

N-(M+1) 

N-2(M+1 ) 

1 

-0.279 

-0.224 

2 

-0.463 

-0.381 

3 

-0.445 

-0.332 

4 

-0.562 

-0.418 

5 

-0.674 

-0.496 

6 

-1.080 

-0 .866 

7 

-1.539 

-1.283 

8 

-2.063 

-1.762 

9 

-2.587 

-2.006 

10 

-2.597 

-2.175 

1 1 

-2.678 

-2.176 

12 

-2.686 

-2.075 

13 

-2.738 

-1.971 

14 

-2.815 

- 1 .764 

15 

-2.839 

- 

16 

-2.824 

- 

17 

-2.833 

- 

18 

-2.810 

- 

19 

-2.799 

■ 

Lower  Bound  Upper  Bound 


Note:  N  is  the  number  of  data  points  in  each  channel  (32  in  this  case) 
M  is  the  number  of  filter  weights 
is  the  error  power 
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i  i“  (  — r  i*"-  >  '  I'l  <1,  J.  .  ( t 

f-  I  '  r 

K  A  y.  I.  j  -‘h  ;  •  IJ'; 

("Ihl-iDA  i  Ll(> 

I  h  (  Ai  h  l  i }  .'I  "i  I  '(..I  I  j  i-J 
i  h  (  ALY- A  — ‘li  IBl  ii'-!  )  I'jt') 

Aiil-  UA  -  I'Y  r-ii 
m  r  A  ;  0>i 
r  X  "Y  y 

I, 1  X  -DY 

lJi.1  1  A;  1=  I ,  |\l 
X  (  I  )  -  A  (  1  '  +AI  IH1  ■  J. 

I  A  Y  1 ,  r  *  ■*  K  *  > ) 
r  .  1- 

U  v'-l'j  ill' 

I II  1  ,  xi,'  I  ■  j  .  !-i 
j  lY  .  (  I  y'H  I  i,Si'-i  I  A  ‘  I  '  .1 

I  I-  (  O'r'  >  I  VO.  Va):  ' 

1  r  (  r  r  “I-  y  ;•  .  U'.';  ■  .  'Vi' 

A!'  Hi  hi.  :!  Hi  'i!  f-  .  , 

(Y  r  H  VY'iIahA 
I i  Ml'lli I'i'l Al’i)-! I  K  1  ■■  i  .  l.H 
J  H--  ■' 

Hr  1'i.jA'i' 

,  o  I  . 

1  i"  i  AH'  I  :  i  '  '  ,•  'il '  ) 

/  I.  -'(IX-:  I  ■  Ai'H'Li '■ 

i’ll  V  i-r-  ,i|  |H|  <  i  .  h  i-:'-,  (  y  >  ,  bAiyV  (  bX  )  .  lii-tHV  ( 
•Ji  ll  H  i  ■  '  i.  •'  Ai..!~  A 

)!  il  r  A  l.iAi  I-  In’-l.'f '.I  r-A— DX.  A  .  h  Y  /  i  li .  r  I'l 
,i  r-  i  riA!  H  A  >  ),  vl- , 

I '  -a:  }■  :  I-  Dr.'  ':■<  !  V  I  '  Yi.r  A  ) 

I  II  H  A- i  ’  t'-UX  +  ( '  +  W 
.1  h  (  AI  H  A  > './AO,  ■//■O,  '/6Ci 
A:  I .  H  A  ■■■  ''  D  V  —  /  + 1' ! )  /  A I r  A 
HI, I  In  '/hO 

AIDA-'  '  +1.1  Y-i  I J  Y  ,  Z  +  UX  +  -'-i-i  iV  / 

A I  r-  f  ’i|  f  'l  I .  r  A  *  A  I  h  H  t  .i  A 

00  *  V'A  l.-j.i'J 

X  (  .1  )  '  X  (  I  )  +  (  I  —  Al.  h  A  )  -iH-i I  ) 

Cril.L  HbSfO) 

J  f-  r --H  X  )  '■■(  n.  VOO,  DU' 

I  ..  ,  D  .0' 

I I. '  I  K, ,  ■  O'  • 

I  I  '  .  '  r  1  -  .  ,  i'* 

1 ‘i-ti  K  m  t  rAi  K  A'  t  T-'i' !  t  H I  • -*-H  (  i  ) 

I  K-  ■'  IJf'M.  f  P  -  /.I 
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i’^pondix  IE  Bur-j's  Coiiplex  Signal.  tlRSA  Conputer  Program 


■OKI  KAN 

00’>:.. 

1.3  (.1(3  :.-;(.i;j4  l'5-l'iAK-8 ,l  rAi.4- 

C: 

NOI  ih  :  ChPBIJiA  C  :>.  1 ,  MI  DA  i>  ;  1  'I'-AI  .IG-80 

>001 

iOl'!  IN) Hr  X 

■  *oo:? 

nKAI .  P  (  3:;' )  ,  Pbh  (32),  PFK  (  32  ) ,  FNk  W  (  32  > ,  GMh  0  <  32  ) 

1003 

rvhM-'-il  U  ( )  j  U'j'L:’  (  )  1  1  j* I i  (  )  <  HAl.f"'  (  )  ■  H  ( i  i'.'  ) 

•004 

KhA!,  Sh  (513) ,  SOflS:. ,  YX(  IdOO) ,  CA!.(  :  OOo) 

1  >00b 

UU!  li-'l  h.X  hi  Pt:  h  j  PLKj  LiLil.P  l.iU'.iil »  P;  Hf’il  ..h'/  H 

•003 

UK  I  1  h  (  6,  O' ) 

;)007 

FL-KNAT  (  1  X,  ■■•.(■■>)■ It-  I.'I-'  HlJli.  LG.  'rNj  bXPAMii,  )J  > 

>003 

KPAiKA,  lOMIUn,  1  G,  '31.!,  PXPAIIU,  DT 

'.<009 

10 

hi  iPnA  l  (215,  3P;.?.  5) 

■:>o  0 

WKJ  I  r.  (6,  15) 

.'01  ;i. 

15 

rOKMAl  (13,  -if- I/-'  I'dIAKI,  rOiUP,  I  h'l. ,  IIU  J  1  ,  12 

(1012 

RF  AD  (  6,  20  )  131  AP  I  ,  13  rop,  NOl  1 1  ,  5?,  I.P 

>013 

20 

Fl'Pi'ir'i  l'  (6)5) 

.'Ol'l 

NK.,:  iP'b,  25)NUM,  LG,  IRTAK  ),  P-.  TOP,  INC,  NGU  :  ,  I.A  'ON 

>015 

J.  1' 

P.:.KliA  rMHl ,  5X,  /1 7,  OR  12  . 

■-<01 3 

CALL.  BFI.l 

(>01  7 

N'  r  'jl.iLKP  > 

■1013 

DPR  I  nr:  R  ,1  LE  2  ( ,  M I ,  LI  I N) »-  X  ) 

0(319 

CALL.  INPUT  (  R  ,  NUM  3N.  RNRL  ■  1.-NR.II ) 

.1020 

WKT  rR.(5,  3(3)NIJM 

002 1 

30 

i"  i .tpi'iA  f  *  F  EP'i  1  A  1  i.*!”'  h’ 1 “  ;  1 ‘rv ) 

0022 

CAI.L.  WRlNUn,  R,  5) 

0023 

I  PI'^8.  -(A  l  AN(  1.  ) 

'3(324 

UAi  1...  RiRI.,.! 

0025 

5(3 

CON  r  ;i  NOP 

0026 

D  R !  •.)  ■■■=>  R 1  -  0  A  ■  1  (  M  C  U  1  -  1  )  *  R.  X  P  A 1 '  1  > 

0027 

RRlJlu  dN  1  (LiRN./H;XPAND--1.  ■ 

(3028 

D- '  60  Kr  >'  1,  PRUI 

(  .  (329 

3R  ( ')••  R  )  ~R  D.iA  1  (  PI —  1  ) ./ I.iR  N— .  5 

0030 

60 

3R  (R'R O-rp]  (^'SRdPR  ) 

0031 

DO  65  I>-/,NLIM 

'  (J032 

RPRIH  I  )-=KR  Al...  (R  (  I  )  ) 

•  .3033 

CNRIJ  (  I  •"•=A  THAI  PRO  1  )  ' 

0034 

6  5 

C(  >M  r  ]  N  JR. 

!  0035 

DO  70  I-'>,NOri 

(3  (3  3'' 6 

?0 

YX  (I  I'^RLOAT  (  I  :■ 

I  003/ 

CAI  1.  NRNPA''l 

0038 

CAl  L  LDPI  Ul  (YX,  PNRl'J,  RIUN,  50,  lOKI,  5(3,  'TOO' 

(3(;K-'9 

CAI  L  R  IGl  (  YX,  CAL,  3-1E.  3) 

(3040 

CAI  1  BhI  1 

i  00^1 

Rip  I  1  E  (  6,  ) 

(30  <(12 

ROKRAT  ( 2.:iv,  I./p  CHi'NMRI.  i  NAVE  -oh'*-) 

0043 

KEAl  Gc',  '->'-'9)NX 

00^  4 

vv  V 

Ro  'sMA  1  (  15  :■ 

CiO^I  5 

oALI.  NKUPAo 

(3(346 

CALL  PPPI  I  'l  (YX,  CNRI--  NUi'L  Ro,  lOiO,  LO:.  ,  /CO) 

(304  / 

CAL  L  R  IGl  (  YX,  (  AL,  ::43,  3) 

0048 

CALL.  DPI  1 

0049 

NKI  I  F  (6,  ■■>  ■0  ) 

>3050 

9  7/' 

RGkNA  1(20/,  1  /-■  CHANNEI  T[  WAL'R:  iOHH /(■  ■' 

0(35 1 

READ  (6,  ■^■-'■-')LL 

(3(352 

K  -0. 

(3(3  5  .'5 

JJU  R>0  I  -  1  ,  MOP 

i  0(354 

80 

p”,-'-i  R  ( 1 )  itiioN.jioa-  ( I ) ) 

/ 
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ORTRAN 

VOo.  J:-: 

00  30;  14  25-MAf<-8- 

P'^'P/ELCA'l  (NI.IM) 

8AVE'=P 

0':>7 

WRITE  (5,  ;E5)p 

'058 

•li* !!' 

1-  oklia  r  ( 2 1 X ,  ••  p  -  ,  r  1  5  7  ) 

'.<59 

LAI -MOM- 1 

"  lAC) 

DO  150  MIL-i,LA;i 

■001 

nnm-nm+i 

'062 

CALL..  Oh‘KL'(h!UMi  PhL#  PhHi  GGUi’^ij  H) 

'0/-. 

P-  <  1 .  -G  ( IvllM  )  -HOuiM.JG  ( >:-  (  Miv  )  )  )  ■'rp 

'064 

HLI.P-P 

065 

:i  00 

IFiMI'IN  Cl  .18  it  Ip)  CO  TO  150 

v-66 

IPdVNM.  L'l  i  81  AP  I  ) CO  TO  15  0 

'067 

lEiMGLK  (NiJN-IS'lAR'i ILIC).  ME.  0)00  Oj  ;  SO 

06:-: 

Lfil-L  oPh (  Nl'j  1  lij  Li  oH  i  KLUl  1  P-  D  \  ) 

069 

ChI  1.  EEI...I 

070 

Col  1...  ):-;l-l...l 

■071 

RL  Ai)  (  6,  'r'9  7  ;■  MX 

'072 

CAI  1.  MEWRAC 

■0/3 

CALI.  KLiPl.i  I  r  (  YX,  8,  NnIJT,  50,  1010,  50,  '00) 

10/4 

WRI  TEO-.,  l::'0)MMM 

■■;)V5 

WRllECj,  120)MtV 

'076 

120 

hL*f‘<r''Al  (L^OXi  *' -i;- -i'r « -if  i\  1 1  Qp  [3. 

'077 

CAI  L  MR  (MM,  C,  15  ' 

•078 

CAi  L  E  IGl  (  YX,  CAL,  1  2,  2) 

»079 

I E  (IP.  EO.  1,  )  CAI .  1 .,  PR  I T  ( MOU  i  ■  D't ,  S ,  Y  X  :• 

^Oc!() 

IE  (EC.  LE.  0)G0  TO  140 

1081 

Rf-AD<6,  9'-'9)MX 

'0:52 

CALI  ALC(8,  MQl.l  1  ) 

.'083 

CALL  MhMRAi-i 

.>0:84 

CAIJ.  WR(NM,G,5) 

085 

CAi  1.  KBPL.01  <  YX,  8,  MiiUT,  1:0,  1010,  50,  ^OO  > 

>0:86 

MRITE  (6,  120)I'.!MM 

>08  7 

CALI  FIGl  ( YX,  CAE,  12,  1  ) 

.'OS  8 

140 

E-HELP 

>089 

150 

i  Oi'l  1  INI  IE 

>090 

Col  L  I3ELI 

>091 

CAI..L  EXI  i 

>(>92 

Ei'E.i 

Kill  IT  INK 

8  CALI.Elj; 

HI-I..L  , 

IMEUT  ,  MR  ,  AT  AM  ,  El  J  lAi  ,  .iMi 

A I  MAG  , 

IMEMRAG,  LBRLOT,  EIGl  ,  Ci'NJG  ,  BPEC 

SPEC  , 

PPM  ,  AI  G  ,  EXIT 

C'P  T  :i  0N:5 

=  ./ijP:  2 

BLOCK 

L.EMG  1  .H 

MAIM. 

8185  (03V762)<^ 

.  $4$$. 

2  ( 000004 ) 

■«--»CCi|-'|HlLEh' - COKI-^H'- 

PHAsh  iJstD  l-Pl-E 

DPI  :LAKAr  JVrS  OOS-SA  14'51? 
EXEi  :n  rAn:|..KS  01263  1411':- 
ASPhliBLY  01  '74  18244 


/ 


t  . 


1 

j 

i 


i 

i 

1 

■< 
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I 

'  '  'F<TKAN  VGA. 


'  >01 
'  )0'/ 

'  >03 
>04 
■  )0b 
>06 
>0/ 

>03 

■i.iOO 

010 

till  b 

>01  2 
■.-13 

■Ml  4  10 

•lib 

■016  20 

■01  / 

■Ml 

'01’?  30 

■■.).<- 0  40 

021 
‘022 


CO  ::i  4  4  2b-l'AK-3l 

SU).-:Ki:iU  l  I NK  I NHUT  ( !< ,  HI  ,  SH.  Y,  L  ) 

CiOl'-VOOlM  INIltX 
COl'ir'LhX  XCi) 

KtOL  V(N ■ ) , / (N- ) 

INOhXtt:! 

RK0;:'(2  ■■■  IIVDF-X  )Y 
IlviUhXtt:2 

6600(2  ■■  I  lv‘))h.X  )  Z 
DO  0=1,  Ml 
Z6! -  .  — Z(  J) 

COM  r  INI  16 
DO  1 0  0  ■  1  ■  F'  1 

X(  I  )=CM6LX  (  Y  (  I  ) ,  Z  (  I  )  )+SIM-;h:M6I  X  (GAUCSC 
CUM  1  II\1IJ6 
Wh'  1  i  fc  (  Ij  I  ..''..I ) 

FOKi'lAT  ( /'.'MX,  CHAIMNi-L  I  ! -X .  '  CHAIMVl- i  T 
DC.!  40  ,i  -1  ,  M  i 
U  6  n  £  ( ,  IC  )/(!),  Z  (  I  ) 

FOKi'iAl  (  "  C  X  6V,  ;i ,  C'X,  K’:-'.  6i) 

CCM I IMOh 

RHI'IJRIV 

FMU 


ROI.1 1  ;i  M£S  i.;:a;  .l.kd: 

CM6LX  ,  GAI.'S3 

OPT  III  NS  ■-  /  OP:  2 

Bl  OCl ,  L.6NG  I  61 

IMPlJI  2/::  (0OiO42>->t- 

.  2  (OOi'OO^  ) 


»^tCONPIL£P  — 
PHAS6 

DECI  APATlVhS 
b  XbCO  I  AfiLhS 
ASSEMBI  .Y 


-  C0R6^ii--:'r 

036 n  6 666 
0062.6'  1 4  Zb'.' 
0('7S6  14f.M2 

012  IS  13S00 


1 

i 


/ 


6A>  ’6 


1  .  i  (.■  > ,  OhL 

I  •  ■ 
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■  RTKAN 

VOX.  lb 

00;  32;  1 1  2':>--l'iAf- 

1)01 

;3I.IEKi;iU  ri,Mr  NK^NNr  O,  IC) 

.)02 

Rb  AL  O (  1  ) 

OOb! 

CGNPI  EX  i3 

>  04 

DO  1  0  I  E  NN 

OO'b 

A-REAI  (G<  I  )  ) 

>06 

)-:-AII'iAi:'i(0(  f  )  ' 

007 

[b  (  IC.  EO.  '.5)WRnb(5.  20) A,  E 

OOP 

IbdC.  EQ,  A)UKr!E(o.  20)A,  B 

009 

20 

bOKNAT  (  2:3 X,  ■"  (  " ,  b  15.  7,  ' ,  ■' ,  b  i:  v  7.  ■  ) 

0  3  () 

10 

CONI INUb 

oil 

RET ORN 

03 

END 

ROUl INE 

S  CAU.ElJ; 

KEC^L  , 

A I  NAG 

Ob'i  ION:-: 

^-./OP;  2 

El  (  II  :K 

LENG'I  H 

WR 

3  50  (000454  )'*■ 

irttCOhl-'ILEK  -  r,  f  I Rh 

PHASh  U3KD  hRhE 

HKCLARAT)  V- S  0070'.* 
b  Xh  CU  I  AOI  HS  00  XP3  1 
ASSbMhU.Y  0109  V  13921 


bORTRAN 

V06.  3  3 

(>0;  32;  31 

'■>00  :i 

bUNCTICN  GAU::;S(  XBAR.  3IGI3A) 

1)002 

R00  1  =  l  ./Pi'RlOr  yATANCl.  )) 

ooo;7' 

10 

corn'  II  lOb 

1)004 

<  -30  ■)■•  (  Rl'':i'U  (  >1  1  —  5) 

0005 

Y  -'’;i  lOT-'-fiXP  ( ->dX./2  :■ 

0006 

YI  RY-4vi:ii:'i'l  -tt-RANd).  3  ) 

0007 

Ib(YIRY  iVI  Y)O0  lO  10 

ooos 

OA;  i:-.9=X#:t:IGNA-i-X)::AR 

0009 

Rb  TURN 

OOlo 

END 

ROUfJIMbS  LAIJ.hD: 

SORT  ,  ATOM  ..  RAN  ,  bXP 


iipriLNP.  -'OP  2 


El  HOP  l,t:;NGiH 

iiAUSP  I'l  \  oOn^:  Z-R )  * 

ijMPILbK - iO.Rb^^* 

PHA'-.E  lISbD  bREb 

DbOl.ARA  I  iVPS  oOA  '7  '  47'. b 
b  X EM  1 1  Ar-::..bS  '>0  ’:-M;  I'I'.yR 

A  E  NO  I,  Y  Oi.0'.:;-:  l:-;9y'j 


ORTKAN 

006  13 

00;  ;32;  56  25-NAH-31  r'A'i'-  1 

001 

SUERGUTINH  EPPIC  (NUI-I,  NN,  PPP  ,  PER,  P  ,  iT,  i33G,  i.'.'Ti,  H  ‘ 

i002 

Oi  I  l4tN3;  I  OiN  H  (  1  1 1  P‘P  P  ( 1 )  i  P'P  P;  )  i  P  ( "  >  <  >  ■'  (  1  1 1  >.’>  iG  (  ■  .4 1  (  1 

'003 

1.  iJl^  Pl  h.X  Lii  PKh  J  PhKi  H»  i-iUGi  GL’I.iHj  P  *  P*i  *  *  ’"M/  P  ).i 

'00/| 

xp.Ri:,"=(0.  ,  0.  1 

N'NMPP-  1 

h,)06 

IP  (N.  NP:.  O'/OG  to  20 

>007 

D'J  10  .J'-^l  ..  MUM 

'003 

PPP  (.J)'"XP.Rn 

'OOS-' 

1 0 

PPR(.J)'=XhRO 

'010 

20 

;3i\!-XPR0 

'01; 

3D  -XPRO 

'012: 

.J  J  -'  iMOPI-N  - 1 

>013 

DO  30 

'014 

Q"P  (..1-t-N-l  )-P'PF(.j;' 

'Olti 

P-P  ( .J  )  -t-^'PR  (  J  ) 

01/:. 

SN  -3,n-K;0!-1.  ig  (  P  )  •«-Q 

01/ 

SD-SD+P'-'^C'.  >N...li  i  (  P  :■  -h2p-CQ1-I..JG  <  Q  > 

'013 

30 

L  UN  r  .LNIJl- 

■01? 

GGiZ  (  NP' ) '--2  •'^SN.'SO 

.'Oi'.O 

I P  (  N.  E'T.  0  )  >■  (  1  )  =iTGG  ( 1  ) 

'0:/i 

IP  >:  N.  EQ.  0  )  GO  TO  ■'■0 

•022 

Dll  40  ..1^1,'4 

)023 

p:'=N-.,J-t 

■o:/4 

H  ■  1  )  ■  ‘L'L'iG  (  'J  )  Ulit''  (  )  '■t'L  Ur'^'Jo  ( • ’Lil.?  (  P. )  ) 

>025 

40 

CCNI'INOP 

'02/-. 

DO  :-'0  ..1=1,0 

.>027 

l:'0 

iTGG  ( .J  ) '  -p-l  ( .J  ) 

•'Yy'S 

GGGl'1  (  MN  )  =500  (  Ni4  ; 

■029 

;i 

>030 

60 

DO  70  ...1-1 ,  ...1..' 

,'03  j 

PP  R  ( .J  )  =  P'P  P  ( . J  )  -t-CC l-i.  .'1-  ( GGG  ( h  IT )  )  PT;.P  ( ..1 )  ->•  P  ( ...'  +  NP'  1  ) 

>032 

PP. P'  ( .J  )  =PP,P  ( .J-t- 1  )  -i-GiZ-G  ( ITN  )  (  PP  R  ( .J+ 1  ;■  ( .J-i- 1  )  1 

/() 

CONT  INPIP" 

:>03.i 

DO  30  I  =  1,I\IN 

'035 

SO 

G  (  1  )  '-GGi  O  I  > 

'033 

Rh  ”  1 JRN 

:>03/ 

P.N)'i 

ROU  1  I  I'lt 

3  CALI. ED: 

CON.JiJ 

OPT  IONS 

=  /TP;  2 

Bl  orK 

LENG'I  H 

BPhC 

(OO/UT'GHi- 

^noJO'IP  II.E- K - 

PHAS 

■  P  0;5P.Li  PRPE 

DPCLAKAT IVES  00622  ITT  56 

PXhCUTABI  hS  00363  14515 

ASShNBl. 

Y  0:120':'  13813 
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DRTRAN 

006.  :13 

00:  34:  02  oo.-MAK 

•001 

SOliROUriNb:;  PRI  I  <|\ICIJ1  .  DT,  3,  YX) 

>002 

RbAI  S(1),YX(1) 

'003 

F  Rh  U 1  1 .  LGAT  ( NOI JT  -  1  )  ,/)jT 

)004 

Y.(  (  1  !'  =— 0.  bO 

■00b 

rti  10  I '"2.0007 

'00/- 

10 

YX ( I )-YX  < I-l )+FRRQl 

•007 

DO  30  J='1.M001 

'008 

WRIlF  (b,  20)  J.  YX(.J).  S(  J) 

'0(1':.' 

20 

FORFIA  I  (  lOX,  13,  bX,  Rib.  /,  bX,  Rib.  /  ■ 

'010 

30 

OONT INUR 

'Oil 

RR  1  UK'IV 

it)  12 

FK!D 

ROU  l  II4t 
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Appendix  II 

On  Multichannel  (Multivariate)  Maximum  Entropy  Spectral  Analysis 


1.  Introduction 

The  univariate  maximum  entropy  spectral  analysis  has  now  been  well  developed 
and  applied  to  many  defense  research  areas  such  as  radar,  sonar  and  geophysics.  There 
has  been  some  work  done  to  extend  the  maximum  entropy  spectral  analysis  to  multivariate 
case.  Whittle  [1]  and  Robinson  [2] [3]  generalized  the  Levinson-Durbin  recursion  to 
the  multivariate  case  by  fitting  both  foward  and  backward  autoregressions  in  a 
stepwise  fashion.  In  this  thesis.  Burg  [A]  has  mentioned  about  the  multichannel  case. 
However  the  computer  programs  for  both  multichannel  and  multivariate  maximum  entropy 
spectral  analysis  were  only  recently  developed  successfully.  Morf  et.al  [5]  developed 
an  algorithm  for  direct  estimation  of  the  normalized  reflection  coefficients  from  the 
observed  data  for  maximum  entropy  spectral  analysis.  They  also  compared  the  spectral 
estimation  with  the  methods  of  Jones  [6],  Nuttall  [7]  and  Strand  [8],  which  are  more 
of  a  direct  extension  of  Burg's  work  to  the  multichannel  (multivariate)  case.  Burg's 
algorithm  does  not  generalize  directly  since  the  forward  and  backward  autoregression 
matrices  are  not  the  same  in  the  multivariate  case,  and  the  forward  and  backward  one- 
step  prediction  error  covariance  matrices  are  different,  although  they  have  the  same 
determinant.  In  this  report,  the  programs  developed  by  Strand  and  Jones  are  applied 
to  real  multichannel  data  and  imagery  data  in  addition  to  a  set  of  test  signals. 

The  merits  of  these  methods  are  closely  examined.  In  spite  of  programming  complexity 
the  multichannel  and  multivariate  maximum  entropy  spectral  analysis  will  have  increased 
application  as  the  real  data  are  almost  always  gathered  in  several  channels.  Data 
from  several  channels  form  a  vector  for  multivariate  study. 

2.  Brief  Mathematical  Analysis 

Let  x^,X2,...,x^  denote  n  zero  mean  vectors  of  dimension  d  each.  The  sample 
estimate  of  covariance  sequence  for  lag  j  is 
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1 

2  n  t+j  t 


(1) 


where  the  prime  denotes  the  transposed  vector.  The  forward  and  backward  predicting 
autoregressions  of  order  p  are  given,  respectively,  as 

^t-k 


=  I  A 
k=l 


-Cb)  = 


'  '  =  I  \ 

^  k»l  ^ 


(2) 


t+k 


where  and  are  d  x  d  matrices,  and  can  be  determined  recursively  [6]  by 


making  use  of  the  estimated  covariance  matrix  in  Eq.  (1) ,  The  recursion  starts  with 

o  o  o 

The  one-step  forward  and  backward  prediction  error  covariance  matrices  are 


(3) 


=  (I  - 

p  P  P  p-1 

s'")  .  (I  -  B'P)A'Pts«>) 

P  P  P  p-1 

The  forward  and  backward  residuals  are,  respectively, 
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The  recursive  equations  are  then  given  by 
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where  U  is  the  sum  of  cross  products  of  forward  and  backward  residuals  at  lag  p, 
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and  V  and  W  are  estimates  of  (n  -  ,  (n  -  respectively, 


n-p 


V  =  I  B, 

t=l 


(p-1) 

t  ' 


(9) 


"  -  T 

t+p  t+p 


(10) 


Although  the  forward  and  backward  autoregression  matrices  and  the  prediction  error 
covariance  matrices  are  different,  the  multivariate  spectra  should  be  identical 
when  calculated  from  the  forward  and  backward  fits  by 


S(f)  =  h[A(f)]“^  [A*(f)] 

or  by 

S(£)  =  h[B(f)  ]"^  [B*(f)  ] 

where 

A(f)  =  I  -  I  A^P) 
k=l 


B(f)  =  I 


-2TTikhf 

e 


h  is  the  sampling  period  and  *  denotes  complex  conjugate  transpose. 

The  above  approach  based  on  the  work  of  Jones  [6]  does  not  guarantee  stability 
and  does  not  generally  produce  a  non-negative  definite  spectrum  as  has  been  pointed 
out  by  Nuttall  [7].  Subsequently,  Nuttall  [7]  and  Strand  [8]  applied  a  weighted 
arithmetic  mean  error  criterion  in  order  to  provide  model  stability  and  to  ensure 
positive  definite  stationary  spectra.  Another  procedure  suggested  by  Morf,  et.al. 
[5]  that  also  meets  the  spectral  requirements  is  to  compute  the  spectrum  from  the 
normalized  reflection  coefficient  matrix  p.  To  obtain  this  matrix,  W  and  V  are 
factored  by  using  Cholesky  decomposition  into  the  product  of  lower  triangular 
matrices  times  their  transposes.  A  new  recursive  procedure  for  and  by 

using  p  in  place  of  Eqs .  (A)  can  then  be  obtained.  Other  recursive  algorithm  has 
been  proposed  [9]  for  the  solution  of  the  normal  equations  for  both  single  and 
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multichannel  data. 


-  37  - 


References  (Appendix  II) 


1.  P,  Whittle,  "On  the  fitting  of  mutlivar iate  autoregressions,  and  the  approximate 
canonical  factorization  of  a  spectral  density  matrix,"  Biometrika,  vol.  50, 

pp.  129-13A,  1963. 

2.  E,  A,  Robinson,  "Mathematical  development  of  discrete  filters  for  the  detection  of 
nuclear  explosions,"  Journal  of  Geophysical  Research,  vol.  68,  pp. 5559-5567 ,  1963. 

3.  R.  A.  Wiggins  and  E.  A.  Robinson,  "Recursive  solution  to  the  multichannel  filtering 
problem,"  Journal  of  Geophysical  Research,  vol.  70,  pp.  1885-1891,  1965. 

4.  J,  P.  Burg,  "Maximum  entropy  spectral  analysis,"  Ph.D.  dissertation,  Dept,  of 
Geophysics,  Stanford  University,  1975. 

5.  M.  Morf,  A.  Vieira,  D.  T.  L.  Lee  an.i  T.  Kailath,  "Recursive  multichannel  maximum 
entropy  spectral  estimation,"  IEEE  Trans,  on  Geoscience  Electronic,  vol.  GE-16,  no. 2 
pp.  85-94,  April  1978. 

6.  R,  n.  Jones,  "Multivariate  autoregression  estimation  using  residuals,"  in 
"Applied  Time  Series  Analysis,"  edited  by  D.  F.  Findley,  Academic  Press,  1978. 

7.  A.  H,  Nuttall,  "Fortran  program  for  multivariate  linear  predictive  spectral 
analysis,  employing  forward  and  backward  averaging,"  NUSC  Technical  Document  5419, 
Naval  Underwater  Systems  Center,  New  London,  Connecticut  1976. 

8.  0,  N.  Strand,  "Multichannel  complex  maximum  entropy  (autoregressive)  spectral 
analysis,"  IEEE  Trans,  on  Automatic  Control,  vol.  AC-22,  no.  4,  pp.  634-640, 

August  1977. 

9.  P,  A.  Tyraskis  and  0.  G.  Jensen,  "Multichannel  autoregressive  data  models," 
submitted  for  publication;  also  presented  at  the  42nd  EAEC  Meeting  in  Taksim- 
Istanbul,  1980. 


-  38  - 


Appe’^.dix  III 


A  BIBLIOGRAPHY  ON  I.iAXIKUM  LNTRCPY  BPLCTHAL 
ANALYSIS  AND  HLLATSD  TSCHNI'^UES 


I.  Introduction 

Recently  there  has  been  strong  research  interest  on  high 
resolution  spectral  analysis  techniques.  This  is  an  important 
area  of  defense  research  because  of  the  numerous  applications  to 
radar,  sonar,  and  geophysical  areas  of  defense  interest.  An 
excellent  publication  is  the  Proceedings  of  the  1976  and  1979  RA^C 
Spectrum  Estimation  ii/orkshop.  Maximum  entropy  spectral  analysis 
is  one  of  a  number  of  high  resolution  spectral  analysis  techniques. 
The  impact  of  the  Burg's  maximum  entropy  spectral  analysis  method 
is  far  more  significant  than  the  technique  itself.  Thus  in  this 
report  we  present  not  only  the  bibliography  of  the  maximum  entropy 
methods  in  one  and  two  spatial  dimensions  but  also  a  number  of 
related  methods  of  high  resolution  spectral  analysis.  One  common 
assumption  with  all  these  methods  is  that  the  data  record  is  short 
and  thus  the  conventional  fast  Fouriter  transform  method  of  spec¬ 
tral  analysis  is  not  suitable.  Probably  because  of  the  snort  length 
record,  the  maximum  entropy  spectral  computation  is  fairly  sensitive 
c  the  presence  of  noise.  In  the  following  sections,  references  are 
arranged  in  the  first  author's  alphabetical  order.  Each  reference 
is  listed  only  once  in  the  report.  Effort  has  been  made  to  provide 
as  complete  list  of  references  as  possible. 
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